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Abstract 

Two new analytical formulations of the acoustic pressure gradient have been devel- 
oped and implemented in the PSU-WOPWOP rotor noise prediction code. The pres- 
sure gradient can be used to solve the boundary condition for scattering problems 
and it is a key aspect to solve acoustic scattering problems. The first formulation 
is derived from the gradient of the Ffowcs Williams-Hawkings (FW-H ) equation. 
This formulation has a form involving the observer time differentiation outside the 
integrals. In the second formulation, the time differentiation is taken inside the in- 
tegrals analytically. This formulation avoids the numerical time differentiation with 
respect to the observer time, which is computationally more efficient. The acoustic 
pressure gradient predicted by these new formulations is validated through compar- 
ison with available exact solutions for a stationary and moving monopole sources. 
The agreement between the predictions and exact solutions is excellent. The for- 
mulations are applied to the rotor noise problems for two model rotors. A purely 
numerical approach is compared with the analytical formulations. The agreement 
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between the analytical formulations and the numerical method is excellent for both 
stationary and moving observer cases. 

Key words: pressure gradient, acoustic scattering, rotor noise, analytical 
formulation, Ffowcs Williams-Hawkings equation, Farassat’s formulation 1A 


1 Introduction 

Acoustic scattering of the noise generated by rotating blades is an area of re- 
search that is not well developed. For example, a helicopter fuselage, a tiltro- 
tor wing, or the duct surrounding a fan, each may substantially modify the 
acoustic signal that arrives at an arbitrary observer location. Such a modifi- 
cation would change both the magnitude and directivity of the acoustic signal 
from what would be observed for an isolated rotor. The effect of a fuselage on 
the noise held generated by a rotating point source was demonstrated by Atalla 
and Glegg [1,2] using a ray-acoustics approach. Laik and Morris [3] showed a 
direct simulation of acoustic scattering by two and three dimensional bodies 
using an extension of the impedance mismatch method. 

Tools exist for predicting fan noise scattering in turbofan engines, but only lim- 
ited work has been done on the acoustic scattering of rotor noise by short ducts 
(i.e. , ducted tail rotors; ducted propellers for compound rotorcraft; ducted fans 
in UAVs, etc.) The various numerical approaches [4-6] to solve the acoustic 
scattering problem use the acoustic velocity on a scattering surface as a bound- 
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ary condition. For example, a rigid surface requires either the satisfaction of 
the impenetrability condition on the surface or zero normal acoustic velocity 
relative to the scattering surface. Most conventional acoustic codes compute 
acoustic pressure at an observer, not the acoustic velocity, but the gradient of 
the acoustic pressure is related to the acoustic velocity through the linearized 
momentum equation. As result, the boundary condition for the stationary 
scattering surface can be written Vp' s ■ n = — Vp' ■ n, where p' is the inci- 
dent acoustic pressure and p' s is the scattered pressure. The calculation of the 
acoustic pressure gradient is, therefore, a key aspect in solving acoustic scatter- 
ing problems. A numerical evaluation of the pressure gradient, which requires 
evaluation of the spatial derivative of acoustic pressure with respect to each 
direction, is the simplest way to calculate the pressure gradient on the surface. 
Nevertheless, it is computationally expensive. Therefore, it is not practical to 
calculate the pressure gradient numerically for a realistic helicopter configu- 
ration, where the scattering computation may require the acoustic pressure 
gradient at thousands or even tens of thousands of collocation points on the 
scattering surface. Furthermore, for the complicated source (rotating blades) 
and scattering surfaces (complete helicopter configuration), it is not easy to 
obtain the pressure gradient numerically. Therefore, it is important to develop 
an analytical formulation for the pressure gradient to enable routine acoustic 
scattering predictions. 


2 Research objective 

The Ffowcs Williams-Hawking (FW-ff ) equation [7] is a powerful tool to solve 
acoustic propagation from arbitrary moving sources such as rotating blades. In 
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this paper, analytic formulations for the pressure gradient are derived starting 
with the FW-H equation for general moving sources and eventually applied to 
rotor noise. 


The analytical formulations have several distinct advantages in terms of nu- 
merical computation. First, no additional input data is needed to predict the 
acoustic pressure gradient beyond what is already required to predict acoustic 
pressure (or at most, numerical differentiation of the input data). Second, the 
retarded time algorithms that will be used have been refined and thoroughly 
tested in various numerical implementations of formulation 1A, which is a 
retarded-time integral representation of the solution of the FW-H equation. 
Finally, by computing the acoustic pressure gradient analytically, rather than 
using a purely numerical approach, significant computation savings (in terms 
of computer run time and memory) and increased robustness are expected. 
Furthermore, the computation of the acoustic pressure from the isolated rotor 
can be computed concurrently with the acoustic pressure gradient. 


The goals of this paper are as follows: 

(1) Develop a computationally efficient analytical formulation for the acoustic 
pressure gradient to provide accurate input data for the boundary con- 
dition for the scattering problems 

(2) Validate the formulation by comparison with available exact solutions 

(3) Apply the formulation to the rotor noise cases 
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3 Acoustic pressure gradient formulations 


The PSU-WOPWOP rotor noise prediction code is used in this work to predict 
the rotor noise (acoustic pressure), as well as the gradient of the acoustic pres- 
sure on the scattering body. The PSU-WOPWOP code is based on Farassat’s 
formulation 1A [11,12], A brief review is given in the next section. 


3.1 Formulation 1A 

Farassat’s formulation 1A [11,12] is an integral representation of the solution 
to the FW-H equation, without the quadrupole source term. It is a retarded- 
time formulation, which can be written as 


where p', p' T , and p' L denote the acoustic pressure, and the monopole source, 
and the dipole source. When the acoustic data surface coincides to the actual 
impenetrable surface, the last two terms become the thickness and loading 
components of the acoustic pressure, respectively. The monopole noise contri- 
bution p’ T can be written as 


p'(x, t) = Pt(x, t) + p' L (x, t), 


( 1 ) 



( 2 ) 



p 0 U n (rM r + c(M r - M 2 )) 
r 2 ( 1 - M r ) 3 


- ret 
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while the dipole noise contribution p' L is written as 



/ 


r 2 (l — M r ) 3 


( 3 ) 


c 


/= 0 


where (x, t) and (y , r) are the observer and source space-time variables, re- 
spectively, r = |x — y| and c is the speed of sound in the undisturbed medium. 
The data surface is described implicitly by the equation /( y, r) = 0, where 
/( y, r) is defined in such a way that V/ = n, which is the unit outward nor- 
mal to the data surface with components rij. The density of the undisturbed 
medium is p 0 . In Eqs. (2) and (3) the subscripts r,n and M imply the dot 
product of the vector with either the unit vector in the radiation direction r, 
outward normal vector n to the surface / = 0, or the surface Mach number 
M, respectively. The dot over a variable indicates source time differentiation. 
The variables Ui and Li are defined by 


where tq are the components of the local flow velocity vector and v t are the 
components of the local blade surface velocity vector and P tJ is the compressive 
stress tensor. Eqs. (4) and (5) are the form used for a permeable surface, which 
is useful if the flow field around the rotor blades becomes transonic — as is the 
case for high-speed-impulsive noise. Eqs. (l)-(3) omit the quadrupole term in 
the FW-H equation, so all significant nonlinear sources should be contained 
within a permeable surface. This enables the inclusion of the contribution of 
those sources without carrying out a volume integration. For an impermeable 
surface, such as the actual blade surface, f/, = and L t = Pijhj. 


Ui = [l- (p/p 0 )]vi + ( pUi/po ), 


( 4 ) 


Lf Pijfij T pUi(u n rq) , 


( 5 ) 


6 


Formulation G1 


Taking the gradient of Eqs. (2) and (3) directly involves complicated algebraic 


manipulations. It is easier to start with the partial differential equation form of 
the FW-H equation and then use the free-space Green’s function to derive the 
new integral formulation. Details of this approach can be found in reference 
13. In this paper, the formulation is revisited with slightly different notation. 

The acoustic pressure gradient can be found by taking the gradient of the 
FW-H monopole and dipole noise terms (neglecting the quadrupole source). 
The gradient of Eq. (1) is 


The next step is to find the acoustic pressure gradient of the monopole and 
dipole noise sources. The governing equation for the monopole noise is 


where 5(f) is the Dirac delta function with support on the data surface / = 0. 
Using the free-space Green’s function 5(g) /Airr, where g = r — t + r/c, the 
monopole component of pressure can be expressed as 


Vp' = Vp' T + Vp' L , 


( 6 ) 



( 7 ) 




— oo —co 


Taking the gradient of Eq. (8) yields 



— OO — OO 


( 9 ) 



— oo — oo 
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where the symbol V x stands for gradient operator with respect to the observer 
variable x. The spatial gradient operator can replaced by a time derivative 
using the relation 


f 5(g) 1 d ( rS(g)\ r 5(g) 


c dt \ r 


( 10 ) 


Combining Eqs. (9) and (10) yields 


4n Vp'r(x,i) = — / / ^Sfl^dydr 


— oo —oo 


t CO 


+ 


r -^S(f)S(9)dydT 


( 11 ) 


— oo — oo 


Using generalized function theory and geometry [14-16] - and following the 
same steps Farassat used in deriving formulation 1A - the gradient of the 
monopole component of the acoustic pressure is found to be 


d [ 

(id r 

mU n 

dt 

icdtJ 
v /= o 

r( 1 - M r ) 


dS 


ret 


[ 

rpoUn 

J 

=0 

_r 2 (l — M r ) _ 


dS 


( 12 ) 


ret 


dE l 

~df 


By recalling that 


and 

it can be easily shown that 
1 


9 \ 1 


1 d \ ] 


dt ^ 

X 

1 - M r dr L ‘ ' ' J 

X 


ret 


f)v C 

— = -(M r f-M), 

or r 


(13) 


(14) 


Ei = ~ J [f E T ] iet dS+ j 
f = o /= o 


(r - M )p 0 U n 
r 2 ( 1 - M r ) 2 


dS , 


ret 


where 


Et — 


Po{U n “h Un) 

_ 1 _ 

' p 0 U n {rM r + c[M r - M 2 )' 

1 

-i 

i 

£ 

to 

1 

ret 

r 2 (l — M r ) 3 


(15) 


(16) 


- ret 


is the combined monopole noise integrand in formulation 1A. Hence, it is al- 
ready available in the noise prediction code. Finally, the monopole component 
of the acoustic pressure gradient can be written 


AnV p' T (-x, t ) 


6 )_ 

Ft 


-c / l f£ ^ S+ / 
/= 0 /= 0 


(r - M)p 0 U n 
r 2 (l - M r ) 2 


- ret 


dS 


(17) 


The observer time derivative of the two integrals can be determined numeri- 
cally. 


The derivation of the gradient of the dipole noise component of the acoustic 
pressure follows the same procedure as used in the monopole noise noise com- 
ponent. The governing equation for the dipole noise is written as 

□Vl = -V ■ M/)], (18) 


thus the dipole component of acoustic pressure is 


t OO 


4? rpi(x,t) = - V 


~d(f)5(g)dydr 


— OO — OO 
t OO 


I j 5(f) L • V x ( 5( fdyd.r\ . 

- rv-S — \ ' ' 


(19) 


— OO — OO 


Using Eq. (10) in the previous integral yields 


H(x,l) = ^ / J y5(f)5(g)dydr+ J J ^5(f)5(g)dydr, (20) 


t OO 


t OO 


— OO — OO 


— OO — OO 


Then if the gradient of the dipole component of acoustic pressure is taken, the 
result is 


4vr Vp' L (x,t) =|J| J J 5(f) L • V ;c dydr 


— OO — OO 
t OO 


r J 


+ 


S(f)L-vJ r ff\ dyd r 


(21) 


— OO — OO 


Note that the observer and the source space-time variables are independent 
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because none of the Dirac delta functions have been used yet in the integration. 
This approach makes it easy to interpret the differential operators. Had the 
integrated results been used, heavy algebraic manipulations would be needed 
and the differential operators would require careful interpretation. 


Using the following relations 

L.v,|UkI' 


= L • V 

L — 2 L r r 


f H9) + ^U( 9 ) 

r I cr 


8{g) - 

cr at 


(22) 


leads to 


L- V, 


r 6(g)' 


= L • V 

L 


\r 
3 L r r 


A <5(9) + hlf (s) 


cr “ 




cr A at 


4 J ^Hfmdydr 

V —OO —OO 


t OO 


+ 


— — \ Lr ^ d(f)S(g)dydr 


— 00 —00 
t 00 


+ 


— — ? Y^d- d(f)S(g)dydT. 


(23) 


(24) 


— 00 — OO 


Again following the procedure used for formulation 1A, Eq. (24) can be rewrit- 
ten as 


4*V J / 1 (x,<)=~{-i 


1 d 


L r r 


/= 0 L 


/=o 


1 — M r dr \yr(l — M r ) / 
L — 3 L r r 


dS 


- ret 


f 

L — 3 L r r 

J 

=0 

r 3 ( 1 - M r ) 


r 2 ( 1 - Af r ) 
dS. 


dS 


ret 


ret 


(25) 


Simplifying Eq. (25) gives the gradient of the dipole noise component of the 
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acoustic pressure 


4tt V p' L (x,t) =^j- / [r£ L ] ret dS + 


/=o 


v /= 0 
L r r — L r M 
r 2 (l - M r ) 2 


[ 

L — L r r 

J 

=0 

_r 2 (l — M r ) _ 




dS + 


ret 


/=0 


ret 

L — 3L r r 

r 3 (l - M r ) 


dS , 


ret 


where Ef, is the combined dipole noise integrand in formulation 1A 


Er =- 


r( 1 - M r ) 2 


+ 


ret 




M 


r 2 (l - M r ) 5 


ret 


L r {rM r + c(M r - M 2 )) 


ret 


(26) 


(27) 


r 2 (l - M r ) 3 

Again, the observer time derivative in Eq. (26) needs to be taken numerically. 


For convenience, Eqs. (17) and (26) are together referred to as formulation 
Gl. This notation parallels that used by Farassat for the monopole and dipole 
formulation 1, which had a observer time derivative outside of the integrals. 
Evaluation of the pressure gradient can now be completed with substantially 
less computational effort than a direct numerical evaluation of the pressure 
gradient. 


Eqs. (17) and (26) have been implemented in the PSU-WOPWOP noise pre- 
diction code to provide the acoustic pressure gradient at an arbitrary observer 
location. The main challenge of this implementation is the calculation of ob- 
server time derivative, d/dt, of the integrals. Care must be taken to ensure 
that the observer position x remains fixed during the calculation of these in- 
tegrals. To simplify the algorithm description, the integrals which must be 
differentiated, surrounded by the braces in Eq. (26), are represented by Q. 
A second-order backward difference algorithm is used to compute the time 
derivative. The general algorithm for the numerical calculation of d/dt is as 
follows: 
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A. Pick T n - n indicates time step and r represents the emission or retarded 
time. 

B. Compute y * (r n ) - each source point is moving, thus at time n, the posi- 
tion of the ?'-th source point is needed. 

C. Save r n , yi ( r n ), velocity, acceleration, etc. for later use as the r n ~ l and 
T n ~ 2 values once n has been incremented. 

D. Compute x {t n ) (based on y l (r n ) and r n ) - if x does not change (i.e. , 
a stationary observer) , then the arrival time t is found explicitly by 
t = r + r/c; if the observer is moving, both the observer position and 
arrival (observer) time must be determined implicitly at the same time. 

E. Calculate Q (y*, r n ; x (t n ) , t™ n ) = Q™ using velocity, acceleration etc. at 


F. Compute t^ n 1 and f 7 n 2 using r n 1 and r" 2 as follows: 


i. t^n 1 = T n 1 + |x (t n ) - Yi (r n 1 )|/c, 

ii. t™n 2 = T n ~ 2 + |x (t n ) - Yi ( r n ~ 2 ) | /c. 

(If x is stationary, the calculation is simpler.) 

G. Calculate Q (y ^ 1 , r n_1 ; x ( t n ) , = Q 7 ^ 1 and Q (y™ -2 , r n_2 ; x (t n ) , 


n — 2 


Q 


H. Calculate 


9Q ^ Qn " — (1 + a YQn + a ( a + 2)Q 


dt a(l + Oi)(t™ — t™ x ) 

for a non-uniform time step. 

I. Interpolate at t*. 

dt 


where a = 


fn +n— 1 

L n L n 


where r denotes source time, t observer time, yi source vector, Xi observer 
vector, c speed of sound, n time index and t* is the specified observer time of 
interest. It is apparent that this procedure is significantly more complicated 
than computing the acoustic pressure. Nevertheless, the additional computa- 
tional effort will be shown to be significantly less than a purely numerical 
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differentiation of the acoustic pressure. 


3.2 Formulation G1A 

The primary drawback of formulation G1 is that numerical time differentiation 
of the integrals is required. If the observer is stationary, then this requirement 
is not a problem because the time history of the integrals can be easily differen- 
tiated numerically. If the observer is moving with respect to the fluid, as in the 
case of a wind-tunnel test, the situation becomes more complicated because 
the formulation requires the observer to be stationary during the evaluation 
of the integrals. Predictions with a moving observer are possible by adjusting 
the observer position at each time in the acoustic-pressure time history; how- 
ever, three evaluation of the integrals are needed to perform a second-order 
difference approximation to the time derivatives at each observer time. These 
extra integral evaluations become unnecessary if the time derivatives are taken 
inside the integrals analytically. 

Although the process of taking the observer time derivatives inside the inte- 
grals and converting them to source time derivatives is not difficult, it is quite 
tedious. The first step is to apply Eq. (13) and then evaluate the source time 
derivatives that results. Some of the key source time derivatives, which are 
the same as Farassat used in the derivation of formulation 1A, are expressed 
as follows: 


dr 


~(M r r-M), 

r 


dr 

dr 


— cM r , 


(28) 

(29) 
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d 

/1\ 1 dr cM r 

(30) 

dr 

\7V r 2 dr r 2 ’ 

d M r 


(31) 

dr 

= -(-M 2 + M 2 ) + M r , 
r 

dL r 

dr 

c 

= L r + ~{M r L r — L m ). 
r 

(32) 


Some new functions are introduced denoted by the following groups of vari- 
ables: 


W = rM r + c(M r - M 2 ), 


r 2 M r - 3crM • M + c(rM r + c(M 2 - M 2 )) 
X\ = — — . 


U (m, n ) = 


r m (l - M r ) n ’ 


V (m, n) = 


dU (m, n ) nrM r + (n — m)cM 2 + mcM r — ncM 2 


dr 


r m +i(i - M r ) n+1 


or 


(33) 

(34) 

(35) 

(36) 


V(m,n ) = nWU{m + 1, n + 1) + c(m — n)M r U(m + 1 ,n). (37) 


These relations will be used in the process of taking the observer time deriv- 
atives inside the integrals in formulation Gl. 


Eqs. (17) and (26) are the starting point for the derivation of formulation 
G1A. Converting the observer time derivative to a source time derivative and 
using the new variables represented by Eqs. (34) and (35), Eqs. (17) and (26) 
become 


47rVp^(x, t) 


1 r r 1 d 
c Jf= o . 1 — M r dr 


t(QU(1,2) + QWU(2,3)) 


'f = o 


d 


1 — M r dr 


(r-M)QE7(2,2) 




(38) 
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4vr Vp' £ (x,t) = - 


i a_ 

c Jf = o L 1 — M r dr 


{f(U r f/(l,2) + (L r -L M )f/(2,2) 


+ 


1 / 

1 d \ 

C Jf = 0 

. 1 — M r dr 1 

1 / 

r i d \ 

C Jf = 0 

[l - M r dr l 


Jf=Ol 

where poU n is defined as Q. 


(L — 3L r r)U(3, 1) 


+ —WU( 2,3)) 
c 

(L — L r r)U (2, 1) 
L r (i — M)U(2, 2) 
dS, 


dS 


- ret 


dS 


- ret 


dS 


- ret 


- ret 


( 39 ) 


Eqs. (38) and (39) can be written in a short hand notation for convenience 

47rVp' r (x, t)=I 1 + I 2 , (40) 

4t rVp' L (x, t) = J 3 + h + h + h, (41) 

where I\ - Iq correspond to each of the integrals in Eqs. (38) and (39). Detailed 
forms of I\ - /g after performing the differentiation of variables with respect 
to the source time are given in appendix A. 


Eqs. (40) and (41), together with the definitions of I\ - Iq, will be referred to as 
formulation G1A and are the main result of this paper. The designation G1A 
is intended to parallel that of Farassat’s formulation 1A, in which the observer 
time derivative is taken analytically inside the monopole and dipole integrals. 
Formulation G1A does not require numerical time differentiation of the inte- 
grals, and, as a retarded-time formulation, is well suited for subsonic source 
motion. Aside from the problem geometry, only the time-dependent input val- 
ues or at most, numerical differentiation of them are required. Furthermore, 
it will be demonstrated with numerical examples that formulation G1A re- 
quires significantly less operations and computer memory than formulation 
Gl. This will be discussed in detail later. The reduction of computational cost 
is important when the formulation is used for the scattering problem. 
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It is worthwhile to point out possible numerical error sources in the developed 
analytical formulations. Numerical errors can be associated with insufficient 
temporal and spatial resolution of the source and a numerical evaluation of 
integrals. An error analysis will be carried out in the following section. Other 
than these error sources, no other numerical errors associated with wave prop- 
agation, such as dissipation or dispersion errors that are important issues in 
Computational Fluid Dynamics (CFD), are involved in the formulations be- 
cause the solution obtained using the free-space Green’s function is exact at 
the far field under the assumption of linear superposition. 

4 Validation of the analytic formulations of the pressure gradient 

The developed analytic formulations will be validated by comparison with 
exact solutions of the pressure gradient for both stationary and moving mono- 
pole sources. In the case of a moving source, both a stationary source in a 
moving stream and a moving source in a stationary stream will be considered. 

4-1 Validation case 1 : a stationary source 

The first validation case is a stationary point monopole source problem. The 
three-dimensional inhomogeneous wave equation is given 


□V(x,f) = q(x,t). 


(42) 


A source which is concentrated at a point is given 


g(x,t) = Q(t)S(x). 


(43) 
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For a point monopole source with a single frequency a the source is given: 


g(x,£) = Ae^h(x), 


(44) 


where A is a complex constant. 


The three-dimensional Green’s function gives the solution of the acoustic pres- 
sure for a point monopole source located at the origin. 

p , (rt) = w^) 

_ A JZ-r/0 («) 

47rr ’ 

where r = |x — y| and t — r/c is called the retarded time. 


Once the pressure is determined, the particle velocity can be calculated from it 
using the linearized momentum equation. The radial component of the equa- 
tion gives 

du r dp' 

dt dr 

The pressure gradient in the r direction is given by 

dp' A 

dr 47t 


iuj 

1 — ^ } e 

cr r 


0 iuj(t—r/c) 


(47) 



Substituting equation (47) into equation (46) gives 


u r (r,t) 


A 

4:7r(iwp 0 ) 



e iu(t-r/c) 


(48) 


Eqs. (45) and (47) give the exact solutions of acoustic pressure and pressure 
gradient at an observer point. In order to validate the pressure gradient for- 
mulations implemented in PSU-WOPWOP, a permeable data surface surface 
enclosing the point source is used. Eqs. (45) and (48) yield the variables U t 
and Li on the surface. In the prediction, A = 47r is used. 
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The spherical permeable surface used has a radius of 0.5m and the polar 
azimuthal angle 6 and polar angle </> are discretized into 36 and 18 panels, 
respectively. 


Fig. 1 shows a comparison of the acoustic pressure at r = 10m for the exact 
solution and prediction. Source frequency is u = lOrad/s. The predicted result 
agrees very well with the exact solution. The order of the error for the peak 
value is less than 1%. This result demonstrates that the permeable surface 
prediction is correct. 


Fig. 2 shows a comparison of the acoustic pressure gradient for the exact 
solution, the predictions with formulation G1 and G1A and the finite difference 
method. It can be seen the predictions match extremely well with the exact 
solution so that each line can barely be differentiated on this graph. This result 
confirms that the analytical formulations of the pressure gradient can be used 
to accurately compute the pressure gradient for a stationary source case. 


Although we demonstrated that the analytical formulations are very success- 
ful, numerical error associated with the grid resolution of the source needs to 
be addressed. Figure 3 shows the pressure gradient prediction for the station- 
ary case depending on the mesh points on the permeable surface. Three test 
cases were used for the grid convergence test: 36 x 18 mesh points, 18 x 9 mesh 
points, 9x4 mesh points. It was found that the numerical error caused by the 
coarse grid resolution reduced the amplitude of the peaks. 
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4-2 Validation case 2 : wind-tunnel case 


The second validation problem is a “wind-tunnel case”, where the source is 
stationary, but the flow is moving with a constant velocity of U. This test case 
is for the validation of the pressure gradient formulations in a moving source 
case. 


The acoustic velocity potential for a stationary source in a uniform stream is 
written 


r 1 D 2 

-&W 2 



a'lx, f) = o, 


(49) 


where D / Dt is the material (or total) derivative operator and (ft' (it, t ) is the ve- 
locity potential. The source and observer are stationary in a uniformly moving 
stream with a Mach number M = Uq/c. 


The Green’s function for the convective wave equation for a harmonic source 
gives the solution for the velocity potential. It is given by, 

AQ 2 

(ft'(x.,t) = h - exp {ik[R — M(3 2 {x — Xg)]}exp{— cot}, (50) 

47rrt 

where 

R = (3\J f3 2 {x - x s ) 2 + {y - y s ) 2 + {z - z s ) 2 , (51) 

and, (3 = l/\/l — M 2 . The retarded time is given by, 

r* =t- ( M cos <d + \Jl- M 2 sin 2 9), (52) 

where R = JJx — x s ) 2 + (y — y s ) 2 + (z — z s ) 2 and cos 0 = (x — x s )/R. 

The acoustic particle velocity can be obtained by the gradient of the velocity 
potential 

v'(x, t) = V<//(x, t). (53) 
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The acoustic pressure is described by the unsteady Bernoulli equation 




po (iix - 


( 54 ) 


The exact solution for the acoustic pressure gradient with respect to the x, y 
and z coordinates is given by, 


dp' 

dx 

dp' 

dy 

dp' 

dz 


= po\iu 


,d£ 
dx 


(. dcP' 

= p 0 IUJ-— 


U d ^'\ 

0 dx 2 ) ’ 

(55) 

rj d V \ 

0 dxdy ) ’ 

(56) 

d 2 ct>'\ 

0 dxdz j ' 

(57) 


Terms that are necessary for calculating the pressure gradient are provided in 
appendix B. 


The procedure of making the pressure gradient predictions is similar to that 
used for the stationary point monopole source. A spherical permeable surface 
enclosing a point source is created and flow passes by the surface with Mach 
number M. The pressure and velocity evaluated on the surface are passed to 
PSU-WOPWOP and used for the prediction of the pressure gradient. Again, 
A = 47T is used in this problem. 


Fig. 4 shows a comparison of the acoustic pressure at a point observer of 
(100.0,0.0,-5.0) for the prediction and the exact solution for Mach number 
M = 0.5 and M = 0.9 cases. The source frequency is u> = lOrad/s. The 
agreement between the FW-H prediction and the exact solution is excellent 
for both low and high Mach number cases. The order of the error for the peak 
is less than 1% for both cases. 


Fig. 5 shows a comparison of the acoustic pressure gradient for the exact 
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solution, the predictions with formulation G1 and G1A and the finite difference 
method. It can be seen the predictions match extremely well, such that the 
exact solution so that each line can barely be distinguished on this graph. 
In addition to validate the methodology, the result implies that the finite 
difference method, which is a purely numerical method, can be used as a 
baseline for the validation of the pressure gradient formulations when the 
source motion is not linear and the exact solution for the pressure gradient is 
not available. 

Validation case 3 : a uniformly moving source 

Now let us consider a monopole source with a strength of q moving uniformly 
in the x— direction with a velocity U in a stationary fluid. The wave equation 
for the acoustic pressure field generated by this moving point source is of the 
form 

□V(x,t) = ^q{t)5( x ~ Ut ) S (v) S ( z ): ( 58 ) 

where q(t) = Ae lu)t . The linearized unsteady Bernoulli equation is given by 

p\x,t) = -Po^r, (59) 

where <f>’ is the velocity potential. Defining if' = —p 0 (j)', one obtains 

□ 2, 0 / (x, t) = q(t)5(x — Ut)5(y)5(z). (60) 

This equation can be solved in many different ways. Many researchers used a 
linear transformation of coordinates analogous to a Lorentz transformation. 
This is given in reference [17] in detail. This transformation enables the re- 
duction of the problem to that of radiation from a stationary source, but it 
involves complicated mathematical manipulations. In the present work, the 
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solution of Eq. (60) is more easily derived by using Farassat’s formulation 1A 
for the integral solution of the FW-H equation. 


Using the free-space Green’s function and the properties of the 5 function, the 
solution of Eq. (60) becomes 


^'(x,f) 


q(t) 

47rr(l — M r ) 


1 

ret 


where q{r) = Aexp(fo;r). 


(61) 


The pressure can be written as 

d 


P = 


dt 


1 d 


<i{t) 



1 — M r dr 

47rr(l — M r ) 

ret 

v(t) 


+ 

ret 

q(r)(cM r — cM 2 ) 

47rr(l — M r ) 2 

47rr 2 (l — M r ) 3 


(62) 


ret 


Eq. (62) is equivalent to the monopole noise term of Farassat’s formulation 
1A for a moving source with a constant velocity. The particle velocity is given 
by the gradient of the velocity potential 

W(M) 


v'(x,f) = V0'(x,t) = 


Po 


(63) 


where 


W' = 


Aiuq{j)V t 


47rr(l — M r ) 


+ g(r)V 


ret 


A 


47rr(l — M r ) 


(64) 


ret > 


The second part of Eq. (64) needs to be evaluated very carefully, ft should be 
noted that 


A 

47rr(l — M r ) 




- ret > 


A 

47rr(l — M r ) 


( 66 ) 


- ret 


Here the left hand side means that the gradient operator is applied after the 
evaluation of the function at the retarded time. In contrast, the right hand 
side of Eq. (65) is carried out before the retarded time is determined. This 
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gives a wrong answer. To avoid making this mistake, it is useful to preform 
the gradient operation before the retarded-time relation is applied (i.e. , before 
the 5(g) term in the free-space Green’s function is integrated). For this case of 
a uniformly moving source, it is useful to express the denominator of Eq. (65) 
in terms of the observer position, initial source position, and the velocity of 
the source. In this form, one obtains 


r(l — M r ) = c(t — t) — M(x i — Ur) 


( 66 ) 


= V ( X “ ^b ) 2 + (! - M2 ){(y - Vs ) 2 + (z - Zs) 2 }, 
where x,y and z are the observer coordinates and x s ,y s and z s are the initial 
source coordinates. Now the problem is given in only x and t variables and 
the explicit dependence on source space-time is eliminated. Now the second 
part of Eq. (64) becomes 


g(r)V 


1 

4v ry / (x - x s ) 2 + (1 - M 2 ){(y - y s ) 2 + (z - z s ) 2 } 


t 


(67) 


Once the exact solutions for the pressure and particle velocity are determined, 
these data on the permeable surface are used to evaluate formulation G1 and 
G1A in PSU-WOPWOP. The same permeable surface used with the previous 
validation cases is used here and A = —4p 0 ir is used to match the strength of 
the monopole source with that for the wind-tunnel case. 


Fig. 6 shows a comparison of the acoustic pressure at an observer located 
(100.0, 0.0, —5.0) for both the prediction and the exact solution for Mach num- 
ber M = 0.5 and M = 0.9 cases. The agreement between the FW-H prediction 
and the exact solution is excellent for both low and high Mach number cases. 
Again, the order of the error for the peak is less than 1%. 


Fig. 7 shows a comparison of the acoustic pressure gradient for the exact solu- 
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tion, the predictions with formulations G1 and G1A and the finite difference 
method. The derivation of the exact solution for the pressure gradient for a 
moving source case requires heavy algebraic manipulations so the exact solu- 
tion obtained from the validation case 2 (the moving stream case) is used since 
both approaches give identical results. The predictions are again in excellent 
agreement with the exact solution. 

These results demonstrate that the pressure and pressure gradient obtained 
from a stationary source in a moving stream are identical with those predicted 
by a moving source in stationary fluid. The validation for a moving source 
source is complete. 

Fig. 8 shows the instantaneous pressure contour for a source wavelength A = 
2m and M = 0.5. The source is moving in the — x direction. The Doppler effect 
of changing wavelength can be seen in the figure. Figs. 9-11 shows contours of 
the pressure gradient. The contours of pressure gradient tend to lean toward 
the direction of the source motion. 


5 Application of the analytic formulations of the pressure gradient 
to rotor noise 

In this section, two representative calculations are performed to demonstrate 
the capability of the new formulations and to provide some indication of the 
efficiency and robustness of the formulations. The first case considers a model- 
scale UH-1H rotor with untwisted blades operating in a non-lifting hover con- 
dition. This test enables simple and fast calculation for both the pressure and 
pressure gradient. The other test case is for the HART-I model rotor in a 
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forward descent flight, which experiences blade-vortex-interaction (BVI) high 
frequency loading on the blades (although the CFD solution does not fully 
capture the BVI). Measured data is not available for the pressure gradient; 
therefore, the predicted pressure-gradient time histories using formulations 
G1 and G1A must be compared with a purely numerical calculation. The fi- 
nite difference predictions are performed by computing the acoustic pressure 
at several points nearby the observer location and then using a second-order 
central finite difference in each of the three spatial directions. 


5.1 Test case 1 : UH-1H model rotor 

A model scale rotor test, conducted by Boxwell, Yu and Schmitz [18]. in 1978 
and later repeated by Purcell [19] in 1988, has been selected for the validation 
of the present analysis and code. The rotor was a one-seventh scale model of 
a UH-1H main rotor with straight, untwisted blades. The model rotor had 
an NACA 0012 airfoil section. The rotor radius R was 1.045 m with a chord 
of 7.62 cm. The model was run at several high-speed hover conditions with 
low thrust. The high-speed hover condition is not of particular interest for 
the validation of the pressure gradient; therefore, a tip Mach number of 0.6 
is selected for the test case. For the hover noise calculation, an Euler solution 
reported by Baeder, Gallman and Yu [20,21] is used as input data. The Euler 
calculation were performed on a C-H grid; only the lower half of the grid was 
used in the CFD calculations by taking advantage of the symmetry of the 
problem. The Euler calculations required approximately 80 min of CPU time 
on a Cray Y-MP. Details of the Euler calculations can be found in references 20 
and 22. 
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Comparisons of the pressure gradient are made for an in-plane microphone 
located 3.09-R from the rotor hub for a stationary observer. Fig. 12 shows 
the total acoustic pressure and the pressure gradient with respect to the x, 
y, and z directions, respectively. The pressure gradient predicted by the two 
analytical formulations are compared to that obtained by the finite difference 
method. The agreement between the analytic formulations and the finite differ- 
ence method is excellent for all components of the pressure gradient. A closer 
examination reveals that the analytical formulations provide much smoother 
results as compared to the finite difference method. The order of the error for 
the peak is 0.1%. 


5.2 Test Case 2: HART-I model rotor 


The forward-flight capability of the new formulations and code is demon- 
strated for a four-bladed rotor representative of the HART-1 model-scale test. 
This case focuses on unsteady blade loading and forward flight. The OVER- 
FLOW CFD code was used to compute the unsteady flow field around the 
rotor. [23,24] A C-mesh topology was been used for the grid with a total grid 
system of 2.4 million points in the near-body region and 15.0 million points in 
the off-body region - in the coarse grid case. The turbulence model used the 
shear stress transport (SST) [25] k — oj by Menter. The rapid dissipation of 
blade-vortex strength makes the prediction of blade-vortex interactions with 
computational fluid dynamics(CFD) difficult. Although the CFD was not fully 
able to capture the BVI loading on the blades — and hence the peaks of pre- 
dicted noise were considerably underpredicted as shown in the references 23 
and 24 — the comparison of the new analytical formulations for pressure gra- 
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client with the finite difference method is still useful to demonstrate its imple- 
mentation in PSU-WOPWOP. 

As in the UH-1H examples, the finite difference result is compared to that of 
analytic formulations to validate the newly developed formulations. For this 
comparison, the observer is located below the rotor plane at a downstream 
position on the retreating side of the rotor. The observer is in motion with 
the rotor to simulate a wind-tunnel test. Although the absolute magnitude 
of the pressure gradient is unknown, confidence in both the derivation and 
implementation of the new formulations would be gained if all of the different 
methods agree. 

Fig. 13 shows the total acoustic pressure and a comparison of the pressure 
gradient at a moving observer for the HART-1 rotor. The analytical formula- 
tions are in a good agreement with the finite difference method. Upon closer 
inspection (not shown), the finite difference result contains a high frequency 
“jitter” that is thought to be of numerical origin. The analytical formulations 
do not exhibit the same “jitter.” In some other cases with a moving observer 
(not shown) it was found that the acoustic pressure gradient predicted by for- 
mulation G1 was sensitive to the method of computing the observer time and 
position. This has not been studied extensively as formulation G1A does not 
suffer in this regard, and also requires less computational effort. 

Table 1 shows a comparison of computational times for formulation 1A (as a 
reference), formulations G1A and Gl, and the finite difference method. The 
finite difference method requires 7 times as much time as formulation 1A but 
formulation G1A only required 3 times as much computation time as formula- 
tion 1A. Formulation Gl requires approximately 5 times as much computation 
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time as formulation 1A or 60 percent more computation time than formulation 
G1A. This demonstrates the significant computational savings of both of the 
analytical formulations and the superiority of formulation G1A. 


6 Concluding remarks 

In this paper, two analytical formulations for the determination of the acoustic 
pressure gradient have been developed and validated by comparison with avail- 
able exact solutions for both stationary and moving point monopole sources. 
It has been demonstrated that the analytical formulations agree very well with 
the exact solution for three different cases. The fact that all three approaches 
give essentially the same results - although they are quite different in ex- 
pression and implementation - gives confidence that both the derivation and 
implementation have been performed correctly. The formulations are applied 
to rotor blades for both hovering and forward-flight conditions. The analytical 
formulations eliminate numerical oscillations, which are present in the finite 
difference method and result in very smooth predictions. 

It has been found that formulation Gl, which evaluates the observer time 
differentiation of the integrals, is a relatively simple formulation but is some- 
what more difficult to implement in PSU-WOPWOP due to the observer time 
differentiation of the acoustic integrals. Furthermore, in at least one case, it 
was found to be sensitive to the choice of numerical algorithm used to find 
the observer time and location. In contrast, formulation G1A, which takes the 
time derivatives inside the integrals, is a somewhat more complicated formu- 
lation. Nevertheless, it yields improved computational efficiency and perhaps 
robustness by avoiding the numerical time differentiation of the acoustic inte- 


grals. Numerical tests show that formulation G1A is the fastest and the most 
efficient algorithm for computing the acoustic pressure gradient. This is im- 
portant for use in calculation of the acoustic scattering, which may require 
several thousand pressure gradient calculations at the collocation points on 
the scattering body. 


A Formulation G1A for the analytical pressure gradient 


Five final equations of formulation G1A for the analytical pressure gradient 
are given as follows: 


h = — 

c Jf = o L 


r{QU(l, 3) + (3QW + QW)U{ 2, 4) + 3QW 2 U(3 , 5)} 


— cM{QU (2, 3) + QWU (3, 4)} 


(A.l) 


dS, 


ret 


1 2 = 


'f = 0 


(M - r)QU( 2, 3) + (~cM r f + cM + rM)QU(3, 3) 


2(M — r)QWU(3, 4) 


(A.2) 


dS, 


ret 


I'A = 


C * Jf = 0 


r(L r + Lr)U ( 1 , 3 ) 

+ c{ — ML r — ( — L r + L]\/[ + L^i^U (2, 3) 

+ £{3L r W + L r W}U (2, 4) 

+ c 2 {(2L r M r — Lm{ 1 + Afr ))f — ( L r — T^f)M }[/ (3, 3) 
+ c{(L r (M r + 2) - 3L M )Wr - L r WM}U(3,4) 


(A.3) 


3L r W 2 rU(3, 5) 


dS, 


- ret 
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n = 


c Jf = o L 


L — L r v)U{2, 2) 

— c{(3L r M r — Lm)y — L r Wl — M r h}U(3, 2) 


(A.4) 


(L — L r r)WU (3, 3) 


dS, 


ret 


I 5 ~ 


l 

c Jf = 0 


{L r (r — M) — L r M}U (2, 3) 

+ c{r {2L r M r — Lm) — M (M r L r — Lm + L r )}U (3, 3) (A-5) 


2L r (r-M)W(3,4) 


dS, 


- ret 


h = 


'/=0 


[L — 3L r r)U(3, 1) 


dS. 


ret 


(A.6) 


Recall that for an impermeable surface, 


Q Po V n ? 


(A.7) 


L Pij Tlj , 


and for a permeable surface, 


(A.8) 


Q p0^n pip^n ^n) 5 


(A.9) 


L Pijflj ~\~ ( 1 1 rj Vj , ) . (A. 10) 

Also note that a dot on the main variables does not imply differentiation of any 
of the associated vectors implied by the subscripts. Subscripts other than i and 
j are a shorthand for the inner product of the main quantity with the vector 
represented by the subscript. The derivative of acceleration, which is called a 
jerk, and second derivative of normal unit vector are evaluated numerically in 
this work. 
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B Exact solution for the pressure gradient for a moving stream 


case 


Terms that are necessary to evaluate the exact solution for the pressure gra- 
dient for a moving stream case are given as follows: 


d4 >' 


dx 


<t>' 


IdR dR 
-R~dx +lk{ &c 
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where 
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Fig. 1. Acoustic pressure comparison for a harmonic stationary source uj = 10 rad/s; 
exact solution: ; prediction with the spherical permeable surface: □ - 



Fig. 2. Acoustic pressure gradient comparison for a harmonic stationary source 

u> = 10 rad/s; exact solution: ; formulation G1A: □ ; formulation 

Gl: o ; finite difference method: o 
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Fig. 3. Acoustic pressure gradient for a different mesh points used; exact solution: 

; prediction with 36 x 18 points: □ ; prediction with 18 x 9 points: 

o ; prediction with 9x4 points: o 


(a) (b) 




Fig. 4. Acoustic pressure comparison for a harmonic source u = 10 rad/s embedded 

in a uniform stream (a) M= 0.5, (b) M = 0.9; exact solution: ; prediction with 

the spherical permeable surface: □ 
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(b) 




Fig. 5. Acoustic pressure gradient comparison for a harmonic source lo = 10 rad/s 

embedded in a uniform stream (a) M= 0.5, (b) M = 0.9; exact solution: ; 

formulation G1A: □ ; formulation Gl: o ; finite difference 

method: o 


(b) 




Fig. 6. Acoustic pressure comparison for a harmonic moving source u) = 10 rad/s (a) 

M= 0.5, (b) M = 0.9; exact solution: ; prediction with the spherical permeable 

surface: □ 


42 


(b) 




Fig. 7. Acoustic pressure gradient comparison for a moving harmonic source u = 10 

rad/s (a) M= 0.5, (b) M = 0.9; exact solution: ; formulation G1A: □ - 

; formulation Gl: o ; finite difference method: o 



Fig. 8. Contour of the acoustic pressure at an instantaneous time for M = 0.5 and 


A = 2 m. The source is positioned at (0,0,0) and the boundary is — 10 < x < 10, 
0 < y < 20, ^ = -5 
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Fig. 9. Contour of dp' /dx at an instantaneous time for M = 0.5 and A = 2 m. The 


source is positioned at (0,0,0) and the boundary is — 10 < x < 10, 0 < y < 20, 
z = - 5 



Fig. 10. Contour of dp' /dy at an instantaneous time for M = 0.5 and A = 2 nr. 


The source is positioned at (0,0,0) and the boundary is — 10 < x < 10, 0 < y < 20, 
z = -5 
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Fig. 11. Contour of dp'/dz at an instantaneous time for M = 0.5 and A = 2 nr. 
The source is positioned at (0,0,0) and the boundary is — 10 < x < 10, 0 < y < 20, 
z = -5 
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Fig. 12. Acoustic pressure and the components of the acoustic pressure gradient for 
the UH-1H model rotor operating in hover with M t ; p = 0.6. (a) total acoustic pres- 
sure (b) dp' /dx (c) dp' /dy (d) dp' fdz\ finite difference method: ; formulation 

G1A: □ ; formulation Gl: o . 
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(a) (b) 






Fig. 13. Acoustic pressure and the components of the acoustic pressure gradient for 
the HART-I rotor operating in a BVI flight condition, (a) total acoustic pressure 

(b) dp' /dx (c) dp' /dy (d) dp' /dz\ finite difference method: ; formulation G1A: 

□ ; formulation Gl: o . 
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Table 1 

Comparison of computational time for the HART rotor with permeable surface 


Formulation 1A 

Formulation G1A 

Formulation G1 

Finite Difference Method 

11.5 (s) 

31.7 (s) 

49.4 (s) 

79.0 (s) 
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